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Abstract. We review recent advancements in modeling the stellar to substellar transition. 
The revised molecular opacities, solar oxygen abundances and cloud models allow to re- 
produce the photometric and spectroscopic properties of this transition to a degree never 
achieved before, but problems remain in the important M-L transition characteristic of the 
effective temperature range of characterizable exoplanets. We discuss of the validity of these 
classical models. We also present new preliminary global Radiation HydroDynamical M 
dwarfs simulations. 
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1. Introduction 

The spectral transition from very low mass 
stars (VLMs) to the latest type brown dwarfs 
is remarkable by the magnitude of the trans- 
formation of the spectral features over a small 
change in effective temperature. It is character- 
ized by i) the condensation onto seeds of strong 
opacity-bearing molecules such as CaH, TiO 
and VO which govern the entire visual to near- 
infrared part (0.4 - 1.2 fim) of the spectral en- 
ergy distribution (hereafter SED); ii) a "veil- 
ing" by Rayleigh and Mie scattering of sub- 
micron to micron-sized aerosols; iii) a weak- 
ening of the infrared water vapor bands due to 
oxygen condensation and to the greenhouse (or 
blanketing effect) caused by silicate dust in the 
line forming regions; iv) methane and ammo- 
nia band formation in T and Y dwarfs; and fi- 
nally v) water vapor condensation in Y dwarfs 
(T e ff < 500 K). Condensation begins to occur 
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in M dwarfs with T eff < 3000 K. In T dwarfs 
the visual to red part of the SED is domi- 
nated by the wings of the NalD and 0.77 jum 
K I alkali doublets which form ou t to as much 
as 200 A from the line center (Alla rd et al.l 
2007b, a). The SED of those dwarfs is there- 
fore dominated by molecular opacities and res- 
onance atomic transitions under pressure (» 
3 bars) broadening condi tions, leaving no win- 
dow onto the continuum dAllardlll990i 1997). 

With the absence of magnetic breaking du e 
to a neutral atmosphere (Moh antv et al.f 2002). 
brown dwarfs should present differential rota- 
tion and clouds should be distributed in bands 
around their surface much a s is shown for 
Jupiter. Burgasseret al. (2002) have suggested 
that brown dwarfs in the L-T transition are af- 
fected by cloud cover disruption. And indeed, 
several objects show even large-scale photo- 
metric variability (Artigau et al. 2009, Radigan 
et al. 2012) — on the order o f 5% to even 
10-30% in the best studied case. iBuenzli et al] 
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d2012l) find periodic variability both in near-IR 
and mid-IR for a T6.5 brown dwarf in simul- 
taneous observations conducted with HST and 
Spitzer. The phase of the variability varies con- 
siderably between wavelengths, suggesting a 
complex atmospheric structure. Recent large- 
scale surveys of brown dwarf variability with 
Spitzer (PI Metchev) have revealed mid-IR 
variability on order of a few percent in > 50% 
of L and T type brown dwarfs. On the basis 
of these results, variability may be expected 
for young extrasolar planets, which share sim- 
ilar Teff and spectral types. And the ubiquity 
of cloud structures in L3-T8 dwarfs strongly 
suggests that these may persist into the cooler 
(> T8) objects. 

The models developed for VLMs and 
brown dwarfs are a unique tool for the char- 
acterization of imaged exoplanets, if they can 
explain the stellar-substellar transition. And 
global circulation models subjected to cloud 
formation in presence of rotation are neces- 
sary to explain the observed weathering phe- 
nomena. Recent ly , Alla rd et al. (2012a,b) and 
Rajpu rohit et al.l d2012l) have published the 
preliminary results of a new model atmosphere 
grid computed with the PHOENIX atmosphere 
code accounting for cloud formation and mix- 
ing from Ra diation HydrpDyna mical (RHD) 
simulations dFrevtag et al.1 120 lOh . In this pa- 
per, we present the latest evolution in modeling 
their SED and their observed photometric vari- 
ability and present new prospectives for this 
field of research. 

2. M dwarfs 

Because oxygen compounds dominate the 
opacities in the SED of VLMs, their syn- 
thetic spectra and colors respond sensibly to 
the abundance of oxygen assumed. Allard et al. 
(2012a,b) compare models using differen t solar 
abundance values dGrevesse et al.llT993l 2007, 
Asplund et al. 2009) and fo und an improved 
agreem ent with constraints (Casagrand e et al.l 
2008a) for the solar abundances obtained us- 
ing RHD simulations by Asplund et al. In 
this paper, we pre sent models based on the 
Caffa uetail (20 lib solar abundances. Taking 
the Greves se et al.l (fl993) solar abundance val- 



ues as reference, the Grevesse et al. (2007) 
results show a reduced oxygen abundance by 
-39%, while Asplund et al. obtain a reduc- 
tion by -34%. These values poses problem for 
the inter pretation of so l ar aste ro-seismological 
results dBasu & Antial 120081 Antia & Basu 
201 1). More recently. ICaffau et alJ d201 ll) used 
the C05BOLD code to obtain a more conser- 
vative reduction of oxygen by -22%. The lat- 
ter value still allows an acceptable representa- 
tion of the VLMs while preserving the astero- 
seismological solar results. Higher spectral 
resolution and up-to-date opacities also con- 
tributed to improve V LMs models compared 
to previous versions (lHauschildt et al.l 1 19991 
Allard et al. 2001). These models rely on lists 
of molecular transition determined ab initio. 
See the review by Homeier et al. elsewhere in 
this journal for the opacities used in the BT- 
Settl models presented in this paper. 

The comparison of the BT-Se ttl PHOENIX 
models based on the ICaffauet al.ld201ll) solar 
abund ances to the low r e solution spectra and 
to the ICasagrande et al.1 d2008bl) temperature 
scale is shown in Figs. Q] and [2] Fig. Q] shows 
an unprecedented agreement with spectral type 
through the M dwarf spectral sequence of the 
models. One can see in Fig. [2] that the new BT- 
Settl models lie s lightly to the blue of the BT- 
Du stv models bv I Allard et all d2012l) based on 
the lAsplund et al.l d2009h solar abundances, but 
largely to t he red of th e AM ES-Cond/Dusty 
models by lAllard et al.1 d200ll) and the BT- 
NextGen models Allard e al.l (l2012l) based on 
the iGrevesse et al.l (11993 ) solar abundances. 
The even lower oxy gen abundance values of 
IGrevesse et all (120071) ca u se the MARCS mod- 
els bv iGustafsson et ail (120081) to lie to the 
right of the diagr am. The NextGen models by 
lHauschildt etafl dj999) also lie to the right of 
the diagram due, among others, to the missing 
and incomplete molecular opacities. 

The BT-Settl models are computed solv- 
ing the radiative transfer in spherical symmetry 
and the convect ive transfer using the Mixing 
Length Theory (Bohm-V itensd[l958l see also 
Ludwig et al. 2002 for the exact formalism 
used in PHOENIX) using a mixing length as 
derived by the RHD con vection simulations 
dLudwig et alJl2002ll2006i Freytag et al. 2012) 
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Fig. 1. Fig. 1 of the article by Rajpurohit et al. (in prep.). The optical to red SED of M dwarfs from M0 to 
M9 observed with the NTT at a spectral resolution of 10.4 A are compared to the the best fitting (chi-square 
minimization) BT-Settl synthetic spectra (dotted lines), assuming a solar composition according to Caffau et 
al. (201 1). The models displayed have a surface gravity of long=5.0 to 5.5 from top to bottom. The best fit is 
determined by a chi-square minimization technic. The slope of the SED is particularly well reproduced all 
through the M dwarfs spectral sequence. However, some indications of missing opacities (mainly hydrides) 
persist to the blue of t he lat e-type M dwarf cases such as missing opacities in the B' 2 £ + <- X 2 S + system of 
MgH bv lSkorve t al. ( 2003) and the opacities are totally missing for the CaOH band near 5500 A. Note that 
chromospheric emission fils the NalD transitions in the latest-type M dwarfs displayed here, and that the 
M9.5 dwarf has a flatten optical spectrum due to dust scattering. Telluric features near 7600 A have been 
ignored from the chi-square minimization. 



and a radius as determ ined by interior mod- 
els (lBaraffeetal.1119981 2003, Chabrier et al. 
2000) as a function of the atmospheric param- 
eters (r e ff, surface gravity, and composition). 
PHOENIX use the classical approach consist- 
ing in neglecting the magnetic field, convec- 
tive and/or rotational motions and other multi- 
dimensional aspects of the problem, and as- 
suming that the averaged properties of stars 
can be approximated by modeling their prop- 
erties radially (uni-dimensionally) and stati- 
cally. Neglecting motions in modeling the pho- 
tospheres of VLMs, brown dwarfs, and planets 



is acceptable since the convective velocity fluc- 
tuation effects on line broadening are hidden 
by the strong van der Waals broadening and 
the important molecular line overlapping pre- 
vailing in these atmospheres. But this is not the 
case of the impact of the velocity fields on the 
cloud formation and wind processes (see sec- 
tion |5]below). 

3. Global RHD M dwarf simulation 

VLMs and brown dwarfs are fully convec- 
tive, and their convection zone extends into 
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Fig. 2. Estimate d 7" eff a nd metallicity (decreasing from lighter to darker tones) for M dwarfs by 
Casagrande et al. (2008a) on the left, and brown dwarfs by Golimows ki et al. (2004) an d Vrba et al. 
(2004) on the right are compared to the NextG en isochrones for 5 Gyrs iBaraffe et al.l i 19981) using model 
atmos pheres by various autho rs: MARCS by iGustafsson et al.l l|2008l). ATLAS9 by |Castelli & Kuruczl 
1 2004J) . DRIFT-PHOE NIX bvlHeTling et al] t2008!) . UCM bv lTsuij d2002h . Clear/Cloudy bvlBurrows etafl 
(2006), NextGen by lHauschildt et alJ d!999l) . AMES-Cond/Dustv by I Allard et all d200lh . the ~ST- 
Cond/Dusty/NextGen models by Allard et al. (2012) and the current BT-Settl models. Some curves labeled 
in the legend can only be seen in the more extended Fig. [4] 



their a tmosphere u p to optical depths of even 
10" 3 dAUardlll990i 1997). Convection is effi- 
cient in M dwarfs and their atmosphere, ex- 
cept in the case of 1 Myr-old dwarfs which 
dissociate H2, l ittle sensitive to th e choice of 
mixing length dAllard et al.1 Il997l) . A mixing 
length value betwee n a - 1/Hp - 1.8 - 2.2 has 
been determined bv lLudwig et all (120021 bOQfjl) 
depending on surface gravity. These simula- 
tions have been recently extended to the late- 
type brown dwarf re gime using C05BOLD by 
iFrevtag et al.l ([2010). In brown dwarfs the in- 
ternal convection zone retreats progressively to 
deeper layers with decreasing r e ff . 

The magnetic breaking known to operate 
in low mass stars is not or less efficiently op- 



erating in fully convective M dwarfs (< M3) 
and brown dwarfs. Brown dwarfs can there- 
fore rotate with equatorial speeds as large as 60 
km/sec or periods of typically 1 .5 hour. Some 
have even been observed with periods nearly 
as low as their breakup velocity (1 hour). In 
comparison, planets of our solar system have 
a larger rotation period, such as 10 hours for 
Jupiter, and 24 hours for the Earth. This rapid 
rotation (as well as their magnetic field) causes 
a suppression of the interior convection ef- 
ficiency leading to a slowed down contrac- 
tion during their evolution, and to larger radii 
then predicted by cl assical evolution models 
(IChabrieret al.ll2007l) . 
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Fig. 4. Same as Fig.[2]but extending into the brown dwarf regime for an age of 3 Gyrs. The region below 
2900 K is dominated by dust formation. The dust free models occupy the blue part of the diagram and only 
at best explain T dwarf colors, while the Dusty and DRIFT models explain at best L dwarfs, becoming only 
redder with decreasing r cff . The BT-Settl, Cloudy and UCM r clit = 1700 K models describe a complete 
transition to the red in the L dwarf regime before turning to the blue into the T dwarf regime, only the 
BT-Settl extends into the Y dwarfs regime (r eff < 500 K). 



The C05BOLD code solves the coupled 
equations of compressible hydrodynamics and 
non-local radiative energy transport on a carte- 
sian grid with a time-explicit scheme. It can be 
used in a local setup (with constant downward 
gravity) to model small patches of the stel- 
lar surface or in a global "star-in-a-box" setup 
(with central gravitational potential) to model 
entire st ars. The "star-in-a-box" se tup has been 
used by ISteffen & Frevtagl (120071) to compute 
global RHD simulations of the scaled-down 
Sun in the presence of rotation, described by 
Coriolis and centrifugal forces. In the cur- 
rent code version, angular-momentum conser- 
vation and the isotropy of the flows at small 
Mach numbers have been further improved. 
This setup has been used, see Fig. [3] to com- 



pute a sequence of radially scaled-down toy 
models of M dwarfs with various rotation pe- 
riods to determine the change in mixing length 
that corresponds to the suppression of convec- 
tion by rotation (Schaffenberger & Freytag, in 
prep.). It is found that at realistic Coriolis num- 
bers, rotation is indeed able to reduce the con- 
vective flux (predominantly towards the equa- 
tor) to contribute significantly to the increase 
in radius as suggested from analysing classical 
models. 



4. Cloud formation in PHOENIX 

To describe the clo ud formation in PHOENIX 
lAllard et al.l (120031) took the approach of us- 
ing a cloud model drawn from an extensive 
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Fig. 3. Global RHD simulation of an M dwarfs 
with solar composition. Obtained parameters are 
T cff = 3 143.44 K, logg = 4.38, and a radius of 
0.003 R Q . The surface pressure is 68227.6 g/cms 2 , 
the surface density is 4.52146e-07 g/cm 3 , the cen- 
tral temperature is 28889.7 K, the central pres- 
sure is 1.67214e+09 g/cms 2 , the central density is 
6.92703e-4 g/cm. The model covers 10 pressure 
scale heights from the center to the surface. The sim- 
ulations use grey PHOENIX opacities in the outer lay- 
ers merged with OPAL data for the interior layers. 
The rotation period of the simulation is 1 hour and 
the simulation is shown equator-on in the reference 
frame of rotation. One can notice a slight oblateness 
of the model at this exaggerated velocity for an M 
dwarf. 

study o f cloud form ation in planetary atmo- 
spheres dRossowll 19781) which compares layer- 
by-layer the timescales of the main processes 
(mixing, sedimentation, condensation, coales- 
cence and coagulation). The mixin g timescales 
are ta ken from RHD simulations (iFrevtag et al.l 
2010, see section [5] below for details). For the 
Allard et al. (2012a,b) pre-release of the BT- 
Settl models, the cloud model was improved 
by a dynamical determination of the super- 
saturation — the ratio of the saturation va- 
por pressure to the compound local gas pres- 
sure (P g (t)/P sv /). In the current pre-release the 
cloud model is further improved by the imple- 
mentation of a grain-size-dependent forward 
scattering, and by accounting for nucleation 
based on cosmic rays studies (Tanaka 2005). 



This latter change allows the cloud model to 
limit its refractory element depletion and form 
relatively more dust grains in higher atmo- 
spheric layers. 

The cloud model is solved from the inner- 
most to the outermost atmospheric layer, de- 
pleting the gas composition due to sedimenta- 
tion gradually from the bottom to the top of the 
atmosphere. This model obtains both the den- 
sity distribution of grains and the grain average 
size per layer, and includes 55 types of grains. 
We do not assume a seed composition since our 
cloud model modifies the equilibrium chem- 
istry iteratively to the nucleation rate limit at 
each atmospheric layers. The cloud compo- 
sition in the photospheric layers varies with 
spectral type from zirconium oxide (ZrC>2) and 
refractory ceramics (perovskite and corundum; 
CaTiC>3, AI2O3) in late-M dwarfs, to silicates 
(enstatite, forsterite, etc.) in early-L dwarfs, 
and to salts (CsCl, RbCl, NaCl) and ices (H 2 0, 
NH 3 , NH 4 SH) in late-T and Y dwarfs. This 
method is coherent with the observed weak- 
ening and vanishing of TiO and VO molecu- 
lar bands (via CaTiC>3, TiC<2, and VO2 grains) 
from the optical spectra of late M and L dwarfs, 
revealing CrH and FeH bands otherwise hid- 
den by the molecular pseudo-continuum, and 
the resonance doublets of alkali transitions 
which are only condensing onto salt grains in 
late-T dwarfs. We solve the Mie equation for 
spherical grains using complex refraction in- 
dex of materials as a function of wavelength 
compiled by the Astrophysikalisches Institut 
und Universitats-Sternwarte in Jena. 

The results are shown in Fig. |4] where 
the new BT-Settl model atmospheres colors 
are interpol ated onto the publis hed theoretical 
isochrones (Baraffe et al. 2003]). New interior 
and evolution models consistent with the BT- 
Settl model atmospheres are in development. 

5. Cloud Formation in C05BOLD 

Cloud formation is included in the C05BOLD 
simulations of VL Ms and brown dwarfs 
(Frevtaa et ai1 l2010l) by assuming a mass den- 
sity of forsterite dust "monomers" initially set 
to its maximum abundance at solar metallic - 
ity, i.e. assuming no settling. Evaporation is ac- 
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Fig. 5. 3D RHD simulations using C05BOLD jFrevtag et aTll2010h of a small box of - from top to right 
to bottom right - 2600K, 2200K, 2000K, 1800K and 1500K atmospheres of log £=5.0, solar metallicity. 
The models have no rotation, except for the 2200K case shown at the top left using an unrealistically small 
rotational period of 15 minutes to highlight the - otherwise negligible - effects on the surface convection 
cells and the atmopheric wave pattern. 
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counted for by applying a rate determined from 
the forsterite saturation vapor curve. Grain 
growth by condensation and sedimentation are 
ac counted for us ing the formulations described 
by Rossow dl978l) . The grains are assumed to 
have sizes up to 1 fim and the correspond- 
ing geometric cross-section is added to the 
opacity bins on-the-fly. Fig. shows local 3D 
"box-in-a-star" type RHD simulations using 
C05BOLD and the same cloud model, where 
the effects of cloud formation (restricted to 
forsterite formation) and Coriolis effects (in the 
top left case) are investigated. Here, the cen- 
trifugal force has been neglected since over the 
box size the centrifugal force would simply 
translate into a reduced gravity of the simula- 
tion. The Coriolis force is exaggerated to ex- 
plore the maximum possible effect of rotation 
locally at the surface. 



Global simulations of quasi-hydrostatic M, 
L and T dwarfs that can resolve the convec- 
tion and cloud formation processes, must be 
significantly scaled down (« 20 so far) in ra- 
dius compared to the real object. Conclusions 
from these RHD simulations must therefore 
be wisely drawn from both types of local and 
global simulations to learn about the cloud sur- 
face coverage of across the M-L-T transition. 
The preliminary conclusions are that the ro- 
tation does not affect remarkably the convec- 
tive shape and size of convective surface cells, 
the related gravity waves and the cloud for- 
mation process at the surface. We know also 
from several global circulation simulations of 
rotating planets (periods around 3 days for 
tidally locked for Jupiter around solar type 
stars) that strong winds and currents are gener- 
ated on larger scales at the surface. These large 
scale features are believed to be responsible 
for the distribution of clouds on the surface of 
Jupiter. The higher rotation velocities of brown 
dwarfs can therefore similarly and easily cause 
comparable patterns at their surface. However, 
it is not yet clear how patchy cloud cover- 
age and variability can occur. Corresponding 
global simulations are being developed. 



6. Summary and Prospectives 

We showed that M dwarfs can be modeled ad- 
equately using up-to-date opacities and the re- 
vised solar abundances bv lCaffauet all (1201 lh 
which preserve the agreement with the results 
from solar astero-seismology dAntia & Basul 
1201 lh . We also showed that VLMs and brown 
dwarfs allow to study the process of cloud for- 
mation in a relatively simple context where ir- 
radiation and ground interaction effects impor- 
tant for planets are not present. This allows to 
identify the fundamental mechanism of cloud 
formation and the determi nation of the velocit y 
field by RHD simulations ( Fre vtag et al.l2010l) . 

We have compared the behavior of the re- 
cently published model atmospheres from var- 
ious authors across the M-L-T spectral tran- 
sition from M dwarfs through L type and T 
type brown dwarfs and confronted them to con- 
straints. If the onset of dust formation is oc- 
curring below r e ff = 2900 K, the greenhouse 
or line blanketing effects of dust cloud forma- 
tion impact strongly (J - K s < 2.0) the near- 
infrared SED of late-M and early L-type at- 
mospheres for 1300 < r eff < 2500 K. The BT- 
Settl models by are the only models to span the 
entire regime from stars to planetary mass ob- 
jects (300 K< r eff < 70, 000 K). In the M dwarf 
range, the res ults appear to favor the BT-Settl 
based on the Caffa u et al.1 (1201 lh solar abun- 
dances versus MARCS and ATLAS 9 models 
based on other values. In the M dwarfs range, 
the BT-Settl models show an unprecedented fit 
quality, even in the V bandpass despite still 
missing or incomplete opacities (in particular 
for the electronic bands of CaOH, and missing 
opacities in the corresponding MgH B-X sys- 
tem bands), above 2500 K. In the brown dwarf 
(and planetary) regime, on the other hand, the 
unified cloud model bv lTsuiil(l2002h succeeds 
in reproducing the constraints, while the BT- 
Settl models also show a plausible transition. 
No models succeed in reproducing perfectly 
the M-L transition between 2500 and 2000 K 
at this stage. This T e ff range is similar to 
that of young (directly observable by imaging) 
and strong ly irradiated pl a nets (Hot Ju piters). 
However. lBonnefov et ail d2010l |2013|) obtain 
reasonable parameters and fit quality to the col- 
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ors and SINFONI near-infrared integral field 
spectra of giant planets. 

Possible further development of the cloud 
model could lie in an improvement of the treat- 
ment of coagulation and of the grain opacities 
in general (e.g. porosity, non-sphericity, mixed 
grain composition rather than adding the con- 
tribution of pure grain species, etc.). 

In this paper, we have presented the first 
preliminary global RHD simulations of an M 
dwarf in presence of rotation. This simulation 
is precise enough in the interior while able to 
resolve the atmospheric convective cells at the 
surface. The strong molecular and atomic res- 
onance line opacities prevailing in M dwarf 
atmospheres lead to a control of their cool- 
ing evolution by the atmosphere, such that in- 
terior and evolution codes must consider the 
atmosp heric structure as their boundary con- 
dition (Bar affeet al]|l998l) . This requires the 
computation of large and fine (Ar e ff= 100 K; 
Alogg < 0.5 dex) model atmosphere grids 
which will however be difficult to be achieve 
by computationally expensive RHD codes like 
C05BOLD. Implicit hydrodynamical simu- 
lations more adequate to address the inte- 
rior convection con ditions are being developed 
dViallet etal.ll201ll) . But these simulations will 
have to resolve the atmosphere and treat the 
radiative transfer and the surface cooling ad- 
equately. RHD simulations do not cover evo- 
lution timescales and cannot replace classical 
stellar evolution and atmosphere models. In the 
meanwhile, therefore, classical ID static inte- 
rior models are being developed using the new 
BT-Settl model atmospheres as surface bound- 
ary conditions, and will be published shortly. 

Global RHD simulation which account for 
rotation are needed to resolve the cloud surface 
distribution and induced potential variability. 
In principle, it is only a small step to go from 
the global RHD simulations of an M dwarf 
to global simulations of late-type M dwarfs, 
brown dwarfs, and gas giant planets that ac- 
count for dust cloud formation and rotation. 
Such simulations are under development and 
to be expected for the end of 2013. 

M dwarfs are magnetically active with 
possibly an important magnetic spot coverage 
due to their more concentrated magnetic field 



lines (same magnetic field strengths over 
smaller radii) than solar type stars. Magnetic 
spots are due to the local suppression of 
convection by the magnetic field lines emerg- 
ing at these points of the stellar surface. 
In the Sun, spots are showing an important 
suppression of the local convection, causing 
an important cooling (almost half the solar 
r e ff). The suppression of convection efficiency 
by rotation and magnetic field is also known to 
lead to a slowed down contraction during their 
evolution, and to larger radii then predicted by 
classical interior models for M dwarfs with 
strong magnetic field s in close-binary systems 
dChabrier et al.ll2007l) . Global Magneto RHD 
(MRHD) simulations which account for rota- 
tion are needed to identify the mixing length 
to use in classical models to compensate for 
this suppression. The non-magnetic global 
RHD models presented in this paper already 
show that rotation alone can have a significant 
impact on the convective efficiency and the 
stellar radius. MRHD simulations are being 
devel oped also using C05BOLD (Nuttoetal. 
1201 2l Freytag et al. 2012,Wedemeyer-B6hmet 
al. 2012). The C05BOLD MRHD simulations 
have not yet been applied to global M dwarf 
models (to investigate the dynamo mechanism 
and the suppression of convection by magnetic 
fields in the interior, or to simulate magnetic 
spots). Fortunately however, we do not expect 
effects as important as for solar spots since 
M dwarf model atmospheres have proven 
to be i nsensitive to the v alue of the mixing 
length (lAllardetal.1119971) unless of course if 
convection is completely suppressed at depth 
over the stellar disk. The BT-Settl model at- 
mospheres and synthetic spectra can therefore 
be used safely to determine the parameters 
of moderately active VLMs, brown dwarfs, 
and planetary mass objects. The BT-Settl 
synthetic colors and spectra are distributed via 
the PHOENIX web simulatofl 
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